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I . INTRODUCTION 



Computer-to-computer communications characterizes the 
fundamental link for message and data traffic in many modern 
military communications scenarios. In some cases the link is 
one of many in an extensive network, in others it may 
represent a single point-to-point communications link. 
Regardless, the objective is to transmit bits with a 
satisfactorily low rate of errors from one terminal to 
another. 

The medium for transmitting the data may be wire or 
optical fiber, that is a low pass channel, or it may be 
radio frequency or acoustic, that is a bandpass channel. In 
either case, the actual signal used to carry the data bits 
must be a properly modulated analog signal with sufficient 
energy and positioned in the frequency spectrum to propagate 
effectively through the available medium. In RF links the 
signal will be further restricted to an assigned frequency 
slot or channel in the spectrum. For network links 
transmission may be restricted to assigned time slots. 

What is needed is a signal modulation format that is 
readily adaptable to a variety of link environments, does 
not require special purpose MODEMS to translate between the 
digital and analog domains, can emulate most existing signal 
modulation formats as well as generate entirely new formats 
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and whose descriptive language is the language of Digital 
Signal Processing (DSP) . Multi-Frequency Modulation (MFM) is 
such a technique and research into its implementation as a 
practical communications method is the subject of this 
report. 

Classical modulation methods for bandpass channels use 
amplitude and/or phase to carry signal information on a 
carrier wave in the channel band. When the information 
source is a finite alphabet, as with data or quantized 
analog sources such as digitized speech or video, then only 
a finite number of signal states are needed to represent, or 
code the source. (The presence of one or the other of two 
frequencies, FSK, is also used to transmit binary data, but 
this is less efficient than phase shift keying, PSK, [ref 
1 ] •) 

In order to accommodate multiple sources, individual 
channels are multiplexed together either in the frequency 
domain, ( FDM) or time domain, (TDM) or both. With multiple 
users on a network, Time Domain Multiple Access, (TDMA) , or 
Frequency Domain Multiple Access, (FDMA) is used to allow 
each data source the opportunity to transmit on the net. 

Modulation and multiplexing equipment for communications 
links carrying many channels is primarily analog circuitry 
and is very expensive to purchase and maintain. Both 
commercial and military communications system developers 
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desire to make their equipment as digitally oriented as 
possible in order to increase its flexibility and to 
automate its testing [ref 2]. This trend is being followed 
at SATCOM earth stations, at military switching centers, 
NAVCOMSTA's for example, on ships and aboard aircraft. 

It is our view that ultimate digitization is achieved by 
having the host computers absorb the functions of modulation 
and multiplexing, demultiplexing and demodulation. MFM is a 
modulation technique that is ideally suited to this task 
because its basic structure is one of time and frequency 
slots. In MFM the signals are directly encoded and 
modulated, and decoded and demodulated using digital signal 
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processing techniques within the host computers. 

Connections to the analog link are made thru D/A and A/D 
chips on plug-in boards over DMA I/O channels. Over eighty 
percent of current digital communications signals, including 
spread-spectrum signals, time and frequency slotted signals, 
plus many entirely new signal formats, are simply generated 
and demodulated using MFM. 

In MFM, signal "packets'* arbitrarily located in the 
frequency spectrum and in time are created as shown in 
Fig.l. Each packet consists of L bauds of K tones. These LK 
subsignals form an orthogonal signal set. Each subsignal may 
be independently modulated with phase and amplitude 
information. 

From a theoretical standpoint, recent research has shown 
that multitone signals are more efficient than single 
carrier signals of the same total bandwidth in a general 
linear channel, particuarily if there is a deep null in the 
frequency response [ref 3]. From a practical point of view, 
we have suggested that MFQPSK signals are ideally suited 
for use in high data rate acoustic burst communications from 
moving platforms using medium frequency, relatively wideband 
(25% bandwidth) links [ref 4]. The appealing thing about MFM 
signals is the simple way they can be generated, 
demodulated, and equalized using standard digital signal 
processing algorithims which are now available for nearly 
all industry standard digital computers in a variety of high 
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level languages. Real time processing at high continuous 
data rates is still difficult, however improvements in this 
area are inevitable. A variety of array processor cards are 
available now as add-on cards and new designs (the NEXT 
computer for example) are incorporating DSP processors into 
the computer's basic architecture. Thus, although current 
technology limits direct bandpass generation of MFM signals 
to the tens to perhaps several hundereds of kilohertz range, 
generation of wideband HF radio signals ( signals in the 
tens of MHz range) should be quite feasible within the 
decade. Of course, the audio range signals we can generate 
now are useful in acoustic communications, an application 
that is being pursued by us in conjuction with the Naval 
Ocean Systems Center in San Diego. 

II. THEORY OF MULTI -FREQUENCY MODULATION 

Generation and Demodulation of MFM 

The following definitions are used in MFM (Please refer 
to Figure 1) : 

T: Packet length in seconds 

AT: Baud length in seconds 

k^: Baud length in number of samples 

L: Number of bauds per packet 

At: Time between samples in seconds 
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f x =l/At: Sampling or clock frequency for D/A and A/D 
conversion in Hz. 

Af=l/AT: Frequency spacing (minimum) between MFM tones. 
K: Number of MFM tones. 



Since At=AT/k x , the sample frequency f^kxAf. 
Consequently, there are a maximum of Y~J2 tones spaced at 
intervals of Af Hz between dc and less than iJ2, the 
Nyquist frequency, that can carry amplitude and phase 
information during each baud. Some, or many, of the tones 
may not be used (or equivalently have amplitudes of zero) 
during any or all bauds of the packet. For example, to 
generate bandpass signals between frequencies f x and f 2 , 
only tones between harmonics k^fi/Af and k 2 =f 2 /Af will be 
allowed non-zero amplitudes. Here, the maximum number of 
tones will be K=k 2 -k 1 +l and the signal bandwidth will be 
W=K*Af. Note that the time bandwidth product of the entire 
signal packet, TW=LAT*KAf=LK, is equal to LK, the total 
number of symbols that can be sent in the packet. 

Let us now provide a mathematical description of the 
signal packet. Let the 1 th baud of the transmit signal be 
described by; 



xju) 

where, 



1 x u (u) 

k 



x Ik( u ) = A^cos (27rkAfu + (p^.) 



0<u<AT. 



( 1 ) 

( 2 ) 
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Here, u is time referenced to the beginning of the baud. 
Actual real time is t=t 0 +lAT+u where t 0 is the time of 
initiation of the 0 th baud, i.e., the beginning time of the 
packet. Now the discrete time signal corresponding to the 
1 th baud is found by sampling (1) and (2) at the sampling 
intervals At=l/f x and is given by; 

x i (n) = £ x lk (n) (3) 

where, 

(n) = A lk cos(27rkn/k x + 0^) ; 0<n<k x -l. (4) 

Here, n is discrete time referenced to the beginning of 
the baud. Note that, in general, k may take on all integer 
values between 0 and k x /2 . We will refer to k as the 
"harmonic number" of the MFM tone of frequency kAf. It is 
useful to note that a baud interval, i.e., a time AT, 
contains exactly k cycles of tone k. Thus, adjacent tones 
differ by one in the number of cycles they make during a 
baud . 

Consider the k x point Discrete Fourier Transform (DFT) of 
(3). It is given by ; 

Xi(k') = l X Uc (k«) (5) 

k 

where. 
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(k ' ) = ^k x A lk {exp(-j0 Ik ) 6 (k'-k) 



+ exp ( j0U() <5 ( k * — (kjj-k) ) ) ; 0<k'<k x -l. (6) 

Thus, it is apparent that the discrete time signal (3) is 
given by the k* point Inverse Discrete Fourier Transform 
( IDFT) of (5) , namely; 

x L (n) = IDFTtX^k) ] (7) 

where from (6) it is clear that, 



X : (k) = ^^A^exp ( — j 01k ) ; 0<k<k x /2 

Xi(k x -k) = ^k x A u ,exp(j0 lk ) ; 0<k<k x /2 



( 8 ) 



To summarize, the 1 th baud is generated by taking the IDFT 
of a complex valued array of length k x . The first half of 
the array is loaded with the amplitudes and phases of the 
tones to be included in the MFM signal at the corresponding 
harmonic numbers. The second half of the array is loaded 
with the complex conjugate of the values in the first half 
of the array at the image harmonics (the image harmonic of k 
is k x -k) . This data structure is shown for k x =16 in Fig. 2. 
The complex conjugate symmetry about the midpoint insures 
that the IDFT, also an array of k x points, will be real 
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valued. The values in this 



array, x L (n) , are the 
discrete time samples of the 
analog transmit signal. 
Clocking these samples out 
thru a D/A converter (they 
are stored as binary numbers 
in the computer) at f x 
samples per second completes 
the generation of the 1 th 
baud. The total signal 
packet is generated by an L- 
fold repetition of this 
process . 
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Figure 2 Data Structure for IDFT 
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A block diagram of the transmitter is shown in Fig. 3. 

At NPS we have implemented this using an IBM PC/AT as a host 
computer. Samples may be clocked out at any clock rate up to 
250 KHz; this limit is set by the maximum byte transfer rate 
of the AT ' s DMA channel [ref 5], The maximum packet size is 
61440 samples; this limit is set, currently, at one sector 
of the AT's memory. We have been using this system to 
generate bandpass signals in the 10-14 KHz band and 16-20 
KHz band using a clock rate of 61440 samples per second. The 
relevant signal parameters are given in Appendix 1. 

Demodulation of MFM is accomplished by a process inverse 
to its generation. Given the time domain signal at the 
receiver x x (u) on the interval 0<u<AT: The signal is 
sampled at f x samples per second and converted to digital 
format with an A/D converter. The k x real values thus 
obtained are loaded into a k* point complex valued array 
(the imaginary parts are set to zero) . This becomes the 
array x x (n) for the 1 th baud. Its 1^ point DFT yields the 
complex valued array X x (k) containing, in its first half, 
the amplitude and phase modulation information, A^ and 0^, 
of the K harmonics employed in the generation of the 
transmit signal. 

Notice that the demodulation operation is linear 
(ignoring for the moment the non-linear nature of the A/D 
conversion) . Also notice that only K out of k x complex 
values are retained for decoding. The upper half of the 
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array is redundant; it contains the complex conjugates of 
the values at the image frequencies. Also, non-used 
harmonics in the lower half of the array will, in the 
absence of noise, contain zero values and can be discarded. 
When they do contain noise values, their discard is 
equivalent to filtering out the unused portion of the 
spectrum at the ouput of the demodulation stage. 

Later in this report, it will be shown that DFT 
demodulation of the signal is equivalent to correlation, or 
matched filter, reception for each of the K tones. It will 
also be shown that the tones are all orthogonal; in discrete 
time on the interval 0<n<k x -l, and in continuous time on the 
interval 0<u<AT. Of course, signals in different bauds are 
orthogonal, since they do not overlap in (real) time. Thus, 
the K tones of the L bauds form a set of LK orthogonal 
signals so the response of a matched filter, or correlator, 
to any of the tones (regardless of its modulation) other 
than the one to which it is matched is zero. The fact that 
the receiver is linear, and a matched filter assures that it 
is optimal for demodulation of MFM, that is it maximizes 
signal-to-noise ratio, in the presence of additive white 
noise. This is so in spite of unknown link attenuation at 
each of the frequencies. We note that the channel need only 
have constant attenuation over narrow bands corresponding to 
the bandwidths of the individual tones, Af, and not over the 
entire bandwidth W. Thus the demodulation is still optimal 
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for a wideband signal propagating thru a channel with 
substantial variation in gain across the band. A block 
diagram of the MFM receiver is shown in Fig. 4. 




Figure 4 Block diagram of receiver 



Properties of MFM Signals 

In this section, a number of important properties of MFM 
signals will be presented. In almost every case we shall be 
referring to a single baud, say the 1 th baud, so, unless 
required for clarity of presentation, the subscript 1 will 
be dropped. 



a. Orthogonality of subsignals. 

In continuous time; 

'0 ; k*i 



AT 

J x k (u)Xi(u)du = 1 

O 



(9) 



k ^ ( A k ) 2 AT ; k=i 



12 



and in discrete time; 



1 0 ; i 

( 10 ) 

^ ) 2 k^ ; k=i . 

Eqns (9) and (10) are easily established using (2) and (4) . 

b. Autocorrelation function (acf) and power spectral density 
(psd) of MFM signals. 

The circular (or periodic) acf of an MFM signal is 
defined as; 



r x (p) = I x (n) x (n©p) ; 0<p<k x -l (11) 

r 

where © signifies a left circular shift. Notice that a 
circular shift of any of the harmonic components of x(n) is 
still a harmonic component of the same frequency but with a 
different phase shift. Combine this observation with the 
orthogonality property (10) and (11) is easily evaluated to 
be ; 



r x (P) 

with, 

r k (p) 



l r k (P) 

k 

% ( A k) 2 k x cos ( 27rkp/k x ) . 



( 12 ) 



(13) 
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